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BiOCui-yS was recently reported to superconduct at T ~ 5.8 K at y « 0.15. Band structure 
calculations indicate that the stoichiometric BiOCuS is a band insulator. In this paper, I show 
that the hole-doped (whether in the virtual crystal approximation or with actual Cu vacancies) 
BiOCuS is on the verge of a ferromagnetic instability (c/. Pd metal), and therefore a conventional 
superconductivity with T c ~ 6 K is quite unlikely. Presumably, the hole-doped BiOCuS is another 
example of superconductivity mediated by spin-fluctuations. 
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Recently, Ubaldini el al- reported that in the quater- 
nary compound BiOCui-j,S superconductivity at T c = 6 
K can be induced by a small Cu deficiency (y ~ 15%), 
although not by an electron doping (when substituting 
O with F). Although another similar study has not con- 
firmed the superconductivity^, this result has triggered 
some interest in the community^. I will show below that 
this interest is well justified, for this material is very un- 
likely to be a conventional s-wave superconductor, and, if 
superconducting, is probably another example of a spin- 
fluctuation mediated superconductivity. 

BiOCuS crystallizes in the well known tetragonal Zr- 
CuSiAs structure^, recently made famous by some Fe- 
based high-temperature superconductors (although there 
is no commonality whatsoever between the latter and the 
material in question). The lattice parameters reported 
in Ref. Q are 3.8726 and 8.5878 A, in Ref. i 3.8691 and 
8.5602 A, and in Ref. i 3.868 and 8.557 A. The inter- 
nal coordinates of Bi and S, respectively, are reported as 
0.14829 and 0.6710 (Ref. i) and 0.151 and 0.648 (Ref. 
1. 

In my calculations I used the structure reported in Ref. 
4. The variation of crystallographic parameters withing 
these limits does not affect any conclusions of this paper. 

In agreement with the previous first principles 
calculations 3,4 I found 5 that in the scalar-relativistic ap- 
proximation the stoichiometric compound is a band in- 
sulator, with a direct gap underestimated compared to 
the measured optical gap 4 ", as common in the density 
functional calculations. Adding spin-orbit interaction on 
Bi reduces the gap even further, from 0.75 to 0.55 eV, 
by affecting the unoccupied states. The states right be- 
low the gap are formed predominantly by the Cu xz and 
yz orbitals, while the empty states above the gap have 
quasi-free-electron character (the interstitial region, be- 
ing «60% in volume, contributes up to 50% of the density 
of states of these bands). In agreement with, especially, 
the full potential all-electron calculations of Ref. [H, I 
found a flat band right at the top of the valence band 
at the Z point (Fig. [T]), essentially not dispersing up to 
half-way between Z and X 6 . This band results in a small 
but very sharp density of states (DOS) peak right below 
the band gap (Fig. [5]), with the total weight of slightly 
less than 0.3 electron per copper, and a maximum DOS 
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FIG. 1: Band structure of BiOCuS near the Fermi level. Note 
a flat band in the Z — A direction. 



of 4.1 states/eV. f.u. at 130 meV below the gap (corre- 
sponding to 0.10-0.12 hole/Cu). One can see that this 
peak indeed originates from the flat band along Z — A, 
by looking at the Fermi surface for Ep set at 120 meV 
below the band gap, where the peak just starts (approx- 
imately 0.07 hole/Cu doping). This feature appears as 
cross-shaped Fermi surface pockets emerging from Z (Fig. 

n. 

Given that Cu d orbitals have sizeable Hund rule cou- 
pling this result suggests that Cu ions under 10-20% hole 
doping should be magnetic. Indeed, ferromagnetic (anti- 
ferromagnetic calculation, however, were invariably con- 
verging to zero magnetization—) calculations for y — 0.2 
in the virtual crystal approximation (0.2 hole per Cu) 
converge to a stable solution with the magnetic moment 
M k, 0.12 — 0.13 hb per Cu, but with hardly any energy 
gain. Correspondingly, fixed spin moment calculations 
(Fig. @| show that in the GGA approximation the energy, 
within the computational accuracy, does not depend on 
magnetization up to M — y \ib ■ Note that this magneti- 
zation would correspond to a half-metallic state similar 
to that in the Co-doped FeS^ (in fact the physics of the 
formation of a magnetic state is quite similar here). This 
means that the calculated spin susceptibility is essentially 
infinite, a situation extremely close to Pd metal, where 
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FIG. 2: Density of states of BiOCuS near the Fermi level. 
Note a peak right below the gap, and the leading contribution 
of the Cu states below, and the interstitial states above the 
band gap. (color online) 



carrier density electron-phonon coupling cannot be too 
strong. 

This leads us to the conclusion that the superconduc- 
tivity in the hole-doped BiOCui_ a S must be unconven- 
tional, and probably magnetically mediated, although it 
is not clear yet what particular pairing symmetry can be 
generated in this material. 

One can question the validity of the virtual crystal ap- 
proximation, given that holes are introduced not through 
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FIG. 3: (color online) Left panel: The Fermi surface of 
BiOCui-^S at Ef 120 meV below the gap, where a sharp 
peak in DOS starts. Note a cross-shaped pocket emerging 
along the Z — A direction, responsible for this peak. For com- 
parison, in the right panel the Fermi surface at Ef = 150 
meV (y ~ 0.15) is shown, where the peak is already fully de- 
veloped. Both Fermi surfaces are colored according the the 
absolute value of the Fermi velocity using the same scale (blue 
is zero). 

GGA calculations also yield infinite susceptibility^. As 
in Pd, or, for that matter, in Fe pnictides, fluctuations 
beyond the mean-field level should suppress magnetiza- 
tion and reduce the susceptibility to a large, but finite 
value. But, again as in Pd, strong spin fluctuations are 
bound to destroy any conventional s-wave superconduc- 
tivity, particularly given that with this relatively modest 
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FIG. 4: Fixed spin moment calculations for BiOCui_j,S for 
y — 0.2. The self-consistent solution is indicated by the ar- 
row. Note that the energy is essentially independent of the 
magnetic moment up to the half-metallic magnetization of y 
/lib/Cu. The inset shows a linear growth of the total energy 
when magnetization exceeds the half-metal limit. 



substitution, but through Cu vacancies. To verify that, 
I have performed supercell calculations using four for- 
mula units with one vacancy, that is to say, B14O4CU3S4, 
corresponding to y = 0.25. The calculations converged 
to a ferromagnetic state with M = 0.2 \xb per formula 
unit (0.8 hb per supercell), again slightly below the half- 
metallic limit of 0.25 /is, thus confirming the validity of 
the virtual crystal calculations. Interestingly, the magne- 
tization was quite delocalized, with Cu ions carrying 0.1 - 
0.13 fiB, and the rest of the magnetic moment distributed 
among the sulfur ions and in the interstitial region. 
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0.05 meV. 

In Ref. 0, unconventional notations for the high symmetry 
points are used; conventional A (n, 0, n) is called R, and 
conventional R (tv, it, it) is called A. I am using the standard 
notations. 

I have tried only a checkerboard antiferromagnetic arrange- 
ment, since I do not see any a priori reason for any other 
antiferromagnetic ordering (the fact that the crystal struc- 
ture coincides with that of pnictide superconductors is not 
an argument at all that this material would assume the 



same exotic magnetic structure). Of course, this does not 
mean that there are no spin fluctuations with a finite wave 
vector; it only tells us that if such fluctuations occur, they 
occur at a wave vector, different from 7r, tt (in the unfolded 
Brillouin zone notation). 

I. I. Mazin, Appl. Phys. Lett. 77, 3000 (2000). 

P. Larson, I.I. Mazin, D.J. Singh, Phys. Rev. B69, 064429 

(2004). 



